Hydrodynamic description of transport in strongly correlated electron systems 
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We develop a hydrodynamic description of the resistivity and magnetoresistance of an electron 
liquid in a smooth disorder potential. This approach is valid when the electron-electron scattering 
length is sufficiently short. In a broad range of temperatures, the dissipation is dominated by heat 
fluxes in the electron fluid, and the resistivity is inversely proportional to the thermal conductivity, 
k. This is in striking contrast with the Stokes flow, in which the resistance is independent of k 
and proportional to the fluid viscosity. We also identify a new hydrodynamic mechanism of spin 
magnetoresistance. 
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Hydrodynamics accurately describes most liquids at 
length scales long compared to the particle-particle 
mean-free path, £, but it is rarely relevant to the electron 
liquid in solids. A hydrodynamic description is based on 
the existence of slow variables associated with conserved 
quantities, while neither the momentum nor the energy 
of the electron liquid in a solid is conserved; electron- 
impurity (ei) and Umklapp scattering violate momen- 
tum conservation and electron-phonon (ep) scattering de- 
grades both the momentum and energy of the electron 
fluid. Consequently, even in relatively clean systems such 
as the electron gas in a semiconductor, the kinetics are 
typically described by the Boltzmann equation, and the 
conductivity is related to the corresponding momentum 
relaxation lengths, £ e i and £ ep 

However, there are circumstances in which the electron 
fluid in a semiconductor device, especially when it has 
very high mobility and is moderately strongly correlated, 
r g > 1, exhibits a range of temperatures and sample pu- 
rity where the electron-electron mean-free path, £ ee is 
small compared to the length scales over which momen- 
tum conservation is violated, £ ee <C £ e i and £ ee <C £ ep - 
Moreover, under most circumstances, Umklapp scatter- 
ing is negligible. In this regime, the electron fluid attains 
local equilibrium on the length scale £ ee , which is short 
compared to the scales at which the conservation laws 
break down, so the dynamics of the electron fluid can be 
treated hydrodynamically. In this paper, we develop a 
theory of electron transport in the hydrodynamic regime. 

Linear resistance of a solid object to a hydrodynamic 
flow was considered by Stokes long ago In this case 
the resistance is proportional to the first (shear) viscos- 
ity of the liquid rj, and independent of the second vis- 
cosity £ and the thermal conductivity k, regardless of 
the compressibility of the liquid. The latter property can 
be traced to the fact that the Stokes flow is isentropic, 
i.e. the equilibrium entropy density of the fluid outside 
the obstacle is co-ordinate independent. In the case of 
strongly correlated electronic systems this is generally 
not the case. In the equilibrium state in the presence 



of a random potential the entropy per electron, So(r), is 
inhomogeneous, and consequently electron flow cannot 
be isentropic. We show below that in this case the re- 
sistance depends on all the kinetic coefficients: n, f and 
«. Moreover, in the ideal fluid limit, ft, rj — > 0, the re- 
sistivity diverges as 1/ft, in contrast to the well known 
D'Alambcrt's paradox in Stokes flow. We also show that 
in the hydrodynamic regime the system exhibits a strong 
spin-dependent magneto-resistance. 

While in the present paper we will not analyze any ex- 
plicit experimental system, there are reasons to believe 
that our theoretical results may be relevant to existing 
experiments involving the highly correlated electron gas 
in semiconductor heterostructures. Recently, low density 
two dimensional electronic systems with high mobility 
have become available, in which r s = V/ Ep 3> I and the 
conductivity is relatively high even at low temperatures. 
(V is the characteristic energy of Coulomb interaction, 
and Ep is the bare Fermi energy) These systems exhibit 
unusual temperature and magnetic field dependencies of 
the resistance. (See Refs. |3j-|5| for a review.) If r s > 1 
at T ~ Ep the electron-electron mean free path is as 
short as an inter-electron distance so the hydrodynamic 
approach should be applicable so long as the correlation 
length of the scattering potential £, is large compared 
to the spacing between electrons. (There are even indi- 
cations that the hydrodynamic regime can be realized in 
GaAs MOSFET's at temperatures as high as 300X0,0].) 

Having in mind the linear resistivity we use the Stokes 
approximation, which neglects the nonlinear terms in the 
hydrodynamic velocity. In the absence of external mag- 
netic field, setting all time derivatives to zero (the sta- 
tionary case), the Stokes equations for a charged fluid in 
the presence of an external potential are 

= Vj, (la) 
= n^{d k a[ k - d t P) - {eEi + cW ), (lb) 
= Tj • Vs + divQ . (lc) 

Here T, no and so denote respectively the temperature 
and position- dependent particle density and entropy per 



particle in equilibrium state; j is the particle current den- 
sity; the electron charge is — e and E is the homogeneous 
external electric field. The linear in j corrections to equi- 
librium quantities are indicated by the tilde sign: P and 
U are the current-induced pressure and self-consistent 
potential. The latter is related to the nonequilibrium 
particle density, h, by the Poisson equation. The dissi- 
pative heat flux, Q, and the viscous stress tensor, a' ik , 
are given by 

e'ik = V (dkv> + div k - 2 $ik divtj + (dm, (2a) 

Q = -«Vf. (2b) 

Here v = j/no is the hydrodynamic velocity, d the di- 
mensionality of space, and T is the nonequilibrium cor- 
rection to the temperature due to the flow. The hydrody- 
namic description of the resistivity is fully determined by 
Eqs. (HJ), ((2|), augmented by the equation of state and the 
Poisson equation. The latter determine the spatial dis- 
tribution of equilibrium density no(r) and entropy so(r), 
which are presumed known below. 

To compute the resistivity p one should solve Eqs. (fTJ- 
[3]) and find the current response to the external electric 
field. Alternatively, for a given distribution of the current 
and temperature, the resistivity can be determined by 
equating the Joule heat, pe 2 (j) 2 , to the dissipation rate 
of mechanical energy 

pe 2 (j) 2 = ^-^Q-VT + a^^, (3) 

where (...) = j ■ ■ ■ d d r/ j d r denotes averaging over 
space. This method is more convenient for perturbation 
theory and we will use it in this paper. 

Let us begin with the simplest case of flow in one di- 
mensional (ID) wire. The continuity equation (]f a| re- 
quires the current density to be uniform, j = (j) = const, 
while the hydrodynamic velocity is given by v = (j)/no. 
It follows from Eqs. (fTc)) and (|2b[) that the total heat 
flux q = Tjso — kO x T is uniform (independent of x). 
Its value is established from the condition that the av- 
erage temperature gradient vanishes, yielding q = 
[T(so/k)} j. (In other words, the Peltier coefficient is 
n = -Q/ej = -(T/e)(s //c)[{l/«;)]- 1 .) The temper- 
ature gradient is given by d x T = (T/k)Ssq j, where 
Ss = [s (1/k) — (sq/k)]/(1/k). Substituting these ex- 
pressions into Eq. (|3]) we obtain for the resistivity 

This equation assumes smooth variation of the disorder 
but does not assume smallncss of the relative variations 
of so and no and is valid in the case when the thermal 
conductivity and viscosity are also position-dependent. 



The first term in Eq. ((4]) depends on the amplitude 
of spatial variations of so(r) rather than their gradient. 
Therefore, for any given set of fluid parameters, the re- 
sistivity is dominated by thermal conductivity for a suf- 
ficiently smooth potential. 

Another remarkable feature of Eq. (j4]) is that its first 
term is inversely-proportional to the thermal conductiv- 
ity. This implies that for an ideal fluid K — > the re- 
sistivity diverges, p — > oo. A qualitative explanation of 
the phenomenon is as follows. In equilibrium the entropy 
per electron sq(t) depends only on the fluid density and 
becomes inhomogeneous in the presence of an external 
potential. The flow of an ideal fluid is adiabatic and 
preserves the original entropy per particle. As a result 
the density dependence of the pressure will be different 
in different elements of the fluid. An adiabatic displace- 
ment of such a fluid from its equilibrium configuration 
will induce temperature gradients. The density change 
due to thermal expansion will create a restoring force, 
which is proportional to the displacement of the fluid, 
rather than to the gradient of the displacement. Conse- 
quently at k = ID adiabatic flow is impossible in linear 
response. This argument shows that the resistivity also 
diverges at k — > in two dimensions. It is easy to see 
from Eqs. (fTc)) and (|2b[) that adiabatic flow (n = 0) is 
allowed only along contours of constant sq, of which only 
a set of measure zero percolate across the system in 2D. 

Equation (j4|) is consistent with the result of a micro- 
scopic calculation of the resistivity of a ID system of 
weakly interacting electrons, Eq. (68) of Ref. d. 

Let us now consider the resistivity in two dimensions 
(2D). In contrast to the ID flow, in 2D current conserva- 
tion does not uniquely determine the spatial dependence 
of the current density. It only implies that the latter can 
be expressed as j = (j) + zx Vt/j, where z is a unit vector 
perpendicular to the plane of flow, and the function ip 
describes spatial variations of the stream function Q ■ 

The problem simplifies in the regime where the disor- 
der potential is weak so that perturbation theory can be 
applied. In this case the relative fluctuations of equilib- 
rium density and entropy are small and the viscosities 
and thermal conductivity may be assumed to indepen- 
dent of coordinates. Let us assume that the spatial fluc- 
tuations of the current density are much smaller than the 
average |Vi/>| <C \(j)\- In this case the gradients of the 
hydrodynamic velocity and temperature are linear in the 
inhomogencity, 

div k = {j k )di{n^), Vf = -Q)Ss . (5) 

Substituting these expressions into Eq. ([3]) we obtain 
for the resistivity in the 2D case, 

ftD = _L/l (4 , )° + ( , + {) (vi) 2 Y (6) 
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Although this expression looks similar to the ID result, 
Eq. ((4]), this formula is applicable only to weakly inho- 
mogeneous flows. Therefore all kinetic coefficients are 
coordinate-independent and correspond to those of the 
disorder-free state, and Ssq = sq— (sq). The origin of the 
first two terms in Eq. is the same as that in Eq. (QJ. 
The additional factor of 1/2 corresponds to the inverse 
number of dimensions. It arises because only the gra- 
dients of s and no along the flow contribute to the re- 
sistance. The term containing the shear viscosity r\ in 
Eq. arises because in 2D the inhomogeneous part of 
the hydrodynamic velocity contains shear flow. 

It is interesting to compare this expression with the 
classical expressions for the resistance to a flow past a set 
of fixed objects. In this case, rather than the smoothly 
varying disorder we have treated, one considers the hy- 
drodynamic flow in the presence of spatially sharp ob- 
jects, along the surface of which "stick" boundary con- 
ditions (v(r) = 0) are applied on the hydrodynamic ve- 
locity. Under these conditions, the Stokes formula for 
the resistance produced by a set of macroscopic objects 
embedded into a 2D liquid is 



The corresponding resistivity of the sample is 



1 



\\nR 2 Ni 



(7) 



Here JVj is the concentration of objects and R is their ra- 
dius. Clearly, this expression is similar to that obtained 
from the second term in Eq. although our expres- 
sion contains the sum (77 + £), while Eq. ([7]| contains 
only rj. However, more importantly, as in the ID case, 
when the disorder is sufficiently smooth, the resistance 
is dominated by the thermal conductivity contribution, 
described by the first term in Eq. ©. 

The pcrturbative result, Eq. ([6]). applies when the 
flow is nearly homogeneous. If the thermal conductiv- 
ity k decreases a strongly inhomogeneous flow will de- 
velop even if the equilibrium density is almost uniform, 
|5no/no| <C 1. This is obvious in the k — > limit, 
where the linear response flow is possible only along the 
so = const lines. At finite but small thermal conductivity 
the current flows primarily in narrow channels localized 
near the contours of constant sq that percolate across the 
whole sample. We can estimate the width of the channels 
x and hence P2D m this limit by computing the rate of 
energy dissipation in Eq. (|3J) for an assumed value of x 
and then minimizing with respect to x. The first term 
in Eq. (|3J) may be estimated as 7^ws (C/^) 3 an d the 

second as T ^ Ss °) ) (j) 2 (x/£) 2 , where £ denotes the corre- 
lation length of the disorder potential. Minimizing the 
sum gives x ~ ^a -1 / 6 , where 



a 



KT] 



» 1. 



(8) 



P2D 



1 l Tn((6s r) 



n(n } 2 e 



(9) 



Note that this limit is likely relevant at relatively high 
temperatures, when k is small, although not so high that 
the inequality r ee 3> r ep is violated. In the opposite limit 
a<l the flow is nearly homogeneous, and the resistivity 
is described by the pcrturbative result, Eq. ((6]). 

We now consider the generalization needed to explore 
a spin mechanism of magnctorcsistance. We assume that 
the component of spin parallel to an applied field H is 
conserved, so the spin per electron, a, is a new hydro- 
dynamic variable. For simplicity, we ignore the orbital 
effects of the magnetic field. 

The hydrodynamic equations in the present case 
should be supplemented by the conservation law for spin, 
which in the linear in j approximation reads 



j • V(T = - div j 



(10) 



Here <7o(r, if) is the equilibrium value of a and j CT is the 
density of spin current relative to the fluid. The latter 
consists of a spin diffusion current and a spin thermo- 
current induced by the temperature gradients. Con- 
versely, the heat flux Q acquires an additional contribu- 
tion in the presence of nonequilibrium spin density gra- 
dients, as required by the Onsagcr principle. We can 
express the heat and spin currents in terms of the ki- 
netic coefficients and gradients of temperature and "spin 
chemical potential" as 



(11) 



Here and 712 (if) = 721 (if), and the diagonal elements 
can be expressed in terms of the thermal conductivity 
K and spin diffusion coefficient D a as 7n = T 2 k, and 
722 = D a 

The expression for the resistivity in terms of the dissi- 
pation rate, Eq. ([3]), must now be replaced with 
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_ _ j_ ( 7n 712 \ 


' VT/T" 


_ 


T \ 721 722 / 





pe 2 (j) 2 



, dvi Q-VT 



ik 



dx k 



T 



ia ■ V/2 C 



(12) 



In the weakly inhomogeneous regime we have (in vector 
notation) (Q,j CT ) = — j (TSso, Sao), where Soq = 00 — 
((To)- Using Eqs. (|TT]) and (|T2|) we obtain the resistivity 



5a 



n 



2e 2 \ K ' \ da J T 

(13) 

where 7 1 is the matrix inverse to 7. The magetore- 
sistance arises not only from the if-dcpendcncc of the 
kinetic coefficients n(H), C(if) and 7(if), but also from 
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the equilibrium quantities Ssq(H) and 5<tq{H). Spin po- 
larization by the magnetic field decreases 5s$, thereby 
decreasing thermal dissipation. On the other hand, it 
induces spatial inhomogeneity of the convectivc spin cur- 
rent j<7o(r). This generates diffusive spin currents, which 
increase the resistance of the sample. 

The hydrodynamic results are very general - the 
physics of the particular electron fluid involved enters 
only through the magnitude and functional dependences 
of the various kinetic coefficients on the magnetic field 
and temperature. In the case of r s < li theoretical 
calculations of these coefficients are under good con- 
trol. (See, for example, Ref. @.) For T < E F , 
where Fermi liquid theory applies, Ssq ~ (T/EF)Sno/no, 
rj ~ C ~ mEp/HT 2 £ ee , and k ~ Ep/hT, where m 
is the electron mass. The first term in Eq. ^ is 
on the order ^--^r(<5^o/ n o) an d the second term ~ 

^■|2^--f|-(^o/ 71 o)i wnerc £ i s the correlation length of 
the disorder potential. At temperatures close to the 
Fermi energy the thermal conductivity contribution in 
Eq. (|4]) is larger than the viscous one by a factor £ 2 no ^> 
1, so p ~ -^(Snl/nl). For T > E F , the electrons 
form a classical gas, so Ssq ~ Shq/uq, k ~ vrT/e 2 , 
? y £ r~j rnvrT/e 2 , where vt ~ s/Tjm is the thermal 
velocity. The thermal conductivity term in Eq. © de- 
creases with temperature ~ ^(^o/^o) an( i the viscous 
one increases, ~ ^-(fap/no) ^ ^ 2 , where ~ noe 2 /T is 
the inverse Debye screening length. 

In clean systems at T = the electron fluid crystal- 

(c) 

lizes for large enough values of r s > 1. In 2D case, 

the best estimate of r$ ~ 40 has been obtained by nu- 
merical simulations under the assumption that there is a 
direct transition between the crystal and the liquid states 
[Toj . Although it has been shown [3, 0| that, rather than 
a direct transition, there must be a sequence of transi- 
tions involving electronic microemulsion phases, it still 

(c) 

seems likely that there is a broad interval r« > r s 3> 1 
where at T = the system is in the liquid state. In 
this case, there is a more involved hierarchy of crossover 
scales, since E F C Sip < 1^, where Ep is the (prob- 
ably strongly renormalized, E* F < Ep) Fermi energy, 
Sip ~ Epy/r^ is the plasma frequency at wave vector or 
order n 1 / 2 , and V ~ Epr s is the typical interaction en- 
ergy between electrons. As a result, there are four charac- 
teristic temperature intervals: 1) For T < E F the system 



is in the Fermi liquid regime, which behaves as above. 
2) For E F < T < Q p there is a semi-quantum regime. 
While there is no established theory in this regime, a 
conjecture concerning the T- dependences of s(T), k, 77, 
and £ was put forward in Refs. [3l. Illl. Il2|. Surprisingly, 
there also does not appear to be published experimental 
data concerning the transport properties of liquid 3 He in 
the corresponding regime, much less for the strongly cor- 
related electron liquid. 3) For Sl p < T < V the system is 
a highly correlated classical fluid, a classic problem about 
which much is known empirically, but which is still a sub- 
ject of ongoing theoretical debate 4) For V < T, the 
system is again weakly interacting classical plasma. 
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